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Abstract 

We show that the chaos in Kuramoto-Sivashinsky equation occurs through period-doubling cas- 
cade (Feigenbaum scenario), in contrast, the chaos in Nikolaevskii equation occurs through torus- 
doubling bifurcation (Ruelle-Takens-Newhouse scenario). 
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We concern the onset of spatiotemporal chaos exhibited by Nikolaevskii equation 
in one-dimensional space with periodic boundary conditions, 

d^(x, t) = -dl[e - (1 + dlf]^ - (d^f. (1) 

Because this class of spatiotemporal chaos appears in many types of physical systems, study- 
ing Eq. (??) is important (see, e.g., the introduction part of [^]). When we introduce a 
variable v defined as v = 2d x ip, we can rewrite Eq. (JIJ): 

d t v(x, t) = -d 2 x [e-(l + d 2 x f]v - vd x v, (2) 

which we use in the following. This equation has the two parameters, the bifurcation pa- 
rameter e and the system size L. We know already that the spatiotemporal chaos in Eq. (J2J) 
occurs supercritically when e increases in sufficiently large system [3]. This can be confirmed 
from the fact that the amplitudes of turbulent fluctuations vanish at the limit e = Q, 0] - 
In this paper, we concentrate how the chaos occurs when the system size increases from 
sufficiently small size, i.e., we elucidate bifurcation route to chaos in Eq. (J2J) increasing the 
aspect ratio. If characteristic wavelengths of spatiotemporal patterns are comparable to 
the system size, then the number of active modes accommodated by the system is of order 
1, i.e., the systems behave in much the same way as systems of a few degrees of freedom. 
Thus, we can conjecture that the bifurcation route to chaos in Eq. (j2J) when the aspect 
ratio increases is period-doubling cascade (Feigenbaum scenario) or intermittency transition 
Pomeau-Mannville scenario) or breakdown of T 2 -torus (Ruelle-Takens-Newhouse scenario) 

].flQ. 

First, we concern, instead of Eq. ( |2j) . a well-known model exhibiting spatiotemporal chaos, 
Kuramoto-Sivashinsky (KS) equation 

d t v(x, t) = —dl(l + dl)v — vd x v. (3) 

Its trivial solution v — is unstable with respect to perturbations having wavenumbers 
q G (0, 1) and stable with respect to perturbations having the other wavenumbers. Since 
the wavenumbers are, in finite-size space L, restricted to being 2irv / L with any integer is, 
all wavenumbers lies out of the unstable band (0, 1) if L is smaller than Lf s = 2ir, and 
then the uniform steady state v = is stable. When we increase L slightly from L^ s ', a 
few modes enter in the band (0, 1), i.e., they destabilize the uniform steady state, and the 




FIG. 1: Periodic attractor obtained from the KS equation with L = 2.8L^ S . 



system exhibits a dynamical state that must be described as a low-dimensional attractor. 
To visualize the attractor, we chose a subspace Re^ — KeA 2 , onto which we project the 
attractor, where A u (t) are the Fourier coefficients defined as 



v(x,t) = ^A u (t)e l 



i2ttux/L 



(4) 
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This choice of subspace is relatively good to see topology of the attractor. When we increase 
L from L^ s , the trajectory (ReAi(t), KeA2(t)) is periodic at first as shown in Fig. ^ Then, 
the attractor loses its symmetry as shown in Fig. El After that, the period of the trajectory 
is doubled as shown in Fig. 01 Finally, the attractor becomes chaotic as shown in Fig. |U The 
return map on a Ponicare section of this chaotic attractor is unimodal, as shown in Fig. [SJ 

which means that the KS chaos occurs through the period-doubling cascade (Feigenbaum 

n *« 

scenario). In Ref.[9(, Y. Kuramoto shows a similar unimodal map for Eq. (jHJ) with another 

boundary conditions v = d\v = 0. 

Now, we elucidate the onset of chaos in Eq. ((21). This equation also has the trivial 
solution v = 0. Its unstable band is {yjl — y/e, \J 1 + y/e) . Thus, the critical system size is 



= 2n/ \J\ + y/e (, which is equal to 5.4766 • ■ ■ when e = 0.1). In the following, we use e = 
0.1. This value is sufficiently small in order to observe characteristic spatiotemporal chaos 
exhibited by Eq. (J2J) in large systems 0. Using Eq. (jH), we choose a subspace ReAi — ImAi, 
which is relatively good choice to see topology of the attractor. The attractor in slightly 
larger system-size than is shown in Fig. El where the trajectory seems to be quasi- 
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FIG. 2: Symmetry-broken attractor obtained from the KS equation with L = 2.87L 
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FIG. 3: Period-Doubled attractor obtained from the KS equation with L = 2.871L 



KS 



periodic. In fact, the return map on a Ponicare section of this trajectory is a closed orbit as 
shown in Fig. which means that the attractor is T 2 torus. When we still increase L, we 
see the torus doubled first, as shown in Fig. |H1 and then chaos, as shown in Fig. El Thus, 
the Nikolaevskii chaos occurs through the breakdown of T 2 -torus (Ruelle-Takens-Newhouse 
scenario) . 

Finally, we comment a future issue. In large systems, the spatial power spectrum of the 
Nikolaevskii equation with e > O(0. 1) is qualitatively indistinguishable from that of the 
Kuramoto-Sivashinsky equation [5|. According to this fact, does the Nikolaevskii equation 
with e > O(0. 1) exhibit the Feigenbaum scenario as the Kuramoto-Sivashinsky equation ? 
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FIG. 4: Chaotic attractor obtained from the KS equation with L = 2.872L^- S . 




FIG. 5: Return map for the chaotic trajectory shown in Fig.0J That trajectory crosses the Ponicare 
section Re^4i = near Re^2 = 40 and near Re^2 = 50 one after another. r n are defined as Re^2 
at the intersections near Re^4_2 = 50. The KS chaos occurs through the period-doubling cascade 
(Feigenbaum scenario) because this map is unimodal. 



If so, how does the Ruelle-Takens-Newhouse scenario change to the Feigenbaum scenario ? 

D. T. gratefully acknowledges financial support by the Japan Society for the Promotion 
of Science (JSPS). 
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FIG. 6: Quasi-periodic attractor obtained from the Nikolaevskii equation with L = 3.97L^ and 
e = 0.1. In this figure, the width of the orbit undulates slightly. However, we confirmed that this 
undulation disappears when we plot the trajectory for a longer time interval. It means that the 
trajectory has two incommensurate frequencies. 




FIG. 7: Return map obtained from the Nikolaevskii equation with L — 3.98-L c and £ — 0.1. v n 
are defined as \A\\ on the Ponicare section Re^i = Irm4i > 0. The closed orbit is a cross section 
of T 2 torus. 
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FIG. 8: Return map obtained from the Nikolaevskii equation with L = 3.9845L^ and e = 0.1. The 
definition of r n is the same as in Fig. We see T 2 -torus-doubling. 




FIG. 9: Return map obtained from the Nikolaevskii equation with L = 3.985L^ and e = 0.1. The 
definition of r n is the same as in Fig. [7| The Nikolaevskii chaos occurs through the T 2 -torus- 
doubling bifurcation (Ruelle-Takens-Newhouse scenario). 
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